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The Fermi-LAT gamma-ray data revealed giant ’bubbles’ of emission above and below the Galac¬ 
tic plane with an energy spectrum significantly harder than seen from other directions. How the 
bubbles connect to the Galactic plane is unclear. Previous analyses masked the Galactic plane 
because of the large foreground. In this paper we use a novel spectral template fit, which allows 
a simultaneous determination of the foreground and the hard bubble-like emission in any direc¬ 
tion. We find that bubble-like emission is not only found in the halo, but it is strongly present 
in the Galactic plane as well with a morphology close to the spatial distribution of the 1.8 MeV 
gamma-ray line from ^^Al, a radioactive nucleus synthesized in SNRs. In addition, the spectral 
shape of this hard component coincides with the predicted spectrum from cosmic rays trapped 
in sources (SCRs) Hence, we propose that the bubble-like emission in the plane has a hadronic 
origin, which arises from SCRs. The bubbles in the halo have the same energy spectrum, which 
suggests that they are outfiows from the plane with the gamma-rays arising from hadronic inter¬ 
actions of protons trapped in a plasma of advected gas. Evidence for advected gas is provided 
by the ROSAT X-ray observations from hot gas in the bubble region. Alternatively, the protons 
may be accelerated in the shock wave of the outflow, thus having the same spectrum as the SCRs, 
which are accelerated in shock waves as well. 

Towards the Galactic center (GC) we observe the 1-3 GeV excess, but observe this excess in 
addition in all regions where there is strong ^^Al production. This excludes the dark matter 
annihilation interpretation. Instead, we propose that the excess is caused by a deficit of the CRs at 
low rigidities in the regions of dense molecular clouds, which are characterized by higher energy 
losses and stronger stellar winds. The ^^Al line is a tracer of such regions. If we introduce a 
break in the cosmic ray injection spectra in the regions of strong ^^Al production to effectively 
take the deficit of low rigidity CRs into account, the excess in the Galactic plane and towards the 
GC both disappear. We thus conclude that the excess in the Galactic center is largely an artefact 
from the excess in the Galactic plane. The correlation between the two originates simply from the 
lines-of-sight towards the GC crossing the Galactic plane. 
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1. Introduction 

The FERMI-LAT gamma-ray telescope [1] has surveyed the gamma-ray sky at energies be¬ 
tween 100 MeV and 100 GeV or even above with unprecedented precision. The main contributions 
to the gamma-rays are well understood: cosmic rays (CRs) interacting with the gas and the pho¬ 
tons in the Galaxy, which leads to production by nuclear interactions combined with a smaller 
leptonic component: bremsstrahlung (BR) and inverse Compton scattering (IC) of electrons on 
the photons of the interstellar radiation field (ISRF) [2]. The Fermi Bubbles were discovered as 
an excess over the expected background [3]. The origin of the bubbles is unclear and many pro¬ 
posals have been made, ranging from hadronic cosmic rays (CRs) interacting with hot gas in the 
halo [4,5], to star bursts [6], to AGN activity in our Galaxy [7-10], to dark matter annihilation [11]. 

Currently, the morphology of the bubbles has not been examined in the Galactic plane (GP), 
because of the large foreground, although such a determination could be helpful for the interpreta¬ 
tion. A spectral template fit to the Fermi data allows a simultaneous determination of the foreground 
and bubble signature in all sky directions, including the GP. The idea is simple: the energy spectra 
of the main contributions to the gamma-ray sky are known from accelerator experiments. If com¬ 
bined with the energy spectra of the CRs and the ISRF the gamma-ray spectra can be calculated for 
each contribution using one of the public propagation codes, like Dragon [12] or Galprop [13,14]. 
These are convenient tools, since they include the cross sections from accelerator experiments and 
take the changes in the CR spectra from propagation and energy losses into account. If the known 
contributions of gamma-ray production (tt^ decay, IC and BR) describe the data one should be able 
to fit the gamma-ray spectra in each sky direction by a linear combination of them. Additional 
contributions will be apparent from a poor fit. For every spatial direction, defined as a cone in 
longitude and latitude (sometimes called pixel), the data consist of 21 energy bins with only a few 
free parameters, namely the normalization of each template. This leads to a strongly constrained 
fit in each direction. First results have been published elsewhere [15]. 

The paper has been organized as follows. Sect. 2 describes the data and the template fit, which 
allows for a fine spatial resolution, even to the level of molecular clouds (MCs), thought to harbor 
the sources of CRs. Inside these sources^ the density of CRs and the gas density in the shock 
wave are high, so the gamma-ray production from the CRs inside the sources, called "Source CRs" 
(SCRs) in Ref. [16], may be significant. Since these correspond to gamma-rays from unpropagated 
CRs they should reveal a hard spectrum expected for shock wave acceleration by SNRs with a 
spectral index between 2.0 and 2.3. [17] Evidence for these SCRs is presented and its connection 
with the Fermi Bubbles, which, interestingly, have the same spectral shape as the SCRs, will be 
discussed. Since this emission happens both in the bubbles and the GP we will call it generically 
1/E'^’^ emission. 

In Sect. 3 the morphology of the 1 /E'^’^ emission is determined. In Sect. 4 the energy of the 
this emission is compared with the gamma-ray energy generated by the averaged rate of SNRs in 
our Galaxy, which shows that the rate of SNRs is high enough to generate all I/E^-^ emission. 

* Speaker. 

^Mainly supernova remnants (SNRs), but other sources, like pulsars, may contribute as well. Since SNRs are 
expected to be dominant, we refer to sources as SNRs. 
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Figure 1: Spectral templates for decays of (a): mesons; (b) Bremsstrahlung; (c): inverse Compton 

scattering. Each panel has normalized templates superimposed for all considered subcones in the region 
lbl<60°, |/| < 60°; (d): the Fermi Bubble energy spectrum from Ref. [18], as indicated by the data points 
inside the shaded area. The latter indicates the spread in the different analyses. The shape of these data is 
well reproduced by a template (black line), calculated from production by protons with a l/£’^-^ spectrum, 
as expected for SCRs. For comparison, the IC and BR contribution from Galactic CRs (instead of SCRs) 
are shown as well. They do not describe the shape of the data points. 

In Sect. 5 we show that the emission has a morphology in the GP identical to the 

skymap of the 1.809 MeV line of ^^Al. [19,20] This radioactive element, which is synthesized in 
SNRs, is a good tracer of SNRs, since these are expected to be the dominant sources. [21] The 
excellent agreement between the two skymaps in the GP, both towards the GC and along the spiral 
arms, provides convincing evidence that the l/E'^’^ emission in the GP indeed originates from 
the SCRs. The bubbles in the halo are then interpreted as outflows from the GC into the halo 
with the CRs trapped inside the plasma. These CRs keep the same spectrum as in the sources, 
either because the mean free path in the plasma is so short, that the CRs do not escape by energy 
dependent diffusion or they are accelerated by the shock of the outflow in the halo, thus obtaining 
the spectrum of diffusively accelerated CRs again. 

In Sect. 6 we confirm the excess in the GC with our template fit. This excess has been widely 
discussed before. [22-43] The most exciting interpretations have been the contributions from dark 
matter (DM) annihilation and/or unresolved sources, like millisecond pulsars (MSPs). However, 
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with our template fit we observe this excess in every region of the GP with a strong production. 
Hence, we call it Galacric plane excess (GPE) instead of GCE. If there is a GPE this will show up 
in the GC as well, since the lines-of-sight towards the GC cross the GP A possible interpretation 
will be discussed. Sect. 7 summarizes the results. 


2. Analysis 


We have analysed the diffuse gamma-rays in the energy range between 0.1 and 100 GeV using 
the diffuse class of the public P7REP_SOURCE_V15 data collected from August, 2008 till July 
2014 (72 months) by the Eermi Space Telescope [1]. The data were analysed with the recommended 
selections for the diffuse class using the Eermi Science Tools (EST) software [44]. The sky maps 
were binned in longitude and latitude in 0.5x0.5° bins, which could later be combined at will. 
The point sources from the second Eermi point source catalogue [45] have been subtracted using 
the gtsrc routine in the EST. The recommended selection of events allows one to take the events 
from misidentified hadrons into account, which are part of the isotropic component provided by 
the Eermi software. 

The gamma-ray flux is proportional to the product of the CR densities, the "target densities" 
(gas or ISRE) and the cross sections, but a template fit lumps the product of these three factors into 
a single normalization factor for each gamma-ray component k, thus eliminating the need to know 
them individually. 

The total flux in a given direction can be described as a linear combination of the various 
contributions with known energy spectra (templates): 


^tot ^^ “ 1 “ ^ ^ “ 1 “ tl4\^Bubble ^ “ 1 “ isotropic ^5 ( 2 . 1 ) 


where the normalization factors Hi determine the fraction of the total flux for a given contribution. 
The factors Ui, and hence the flux of each contribution, can be found from a x'^ which tries to 
adjust the templates to best describe the data. Since the number of data points in each spectrum for 
a given subcone (21 energy bins) is large compared with the number of free parameters (nt < 5) 
the fit is strongly constrained, thus allowing a determination of the various template contributions 
in each direction. 

As test statistic we use the x^ function defined as: 


N 21 






n{i^k) xtemplate{ij^k))^ 


( 2 . 2 ) 


where the sum is taken over the N subcones in different sky directions /, data(i,j) represents the 
total Eermi flux in direction i for energy bin j, template(i,j,k) with normalization n{i^k) is the 
contribution of template k out of a total of m templates to data(i,j) and o{iJ) is the total error of 
data(i,j), obtained by adding the statistical and systematic errors in quadrature. The recommended 
systematic errors in the Eermi Software on the total flux are 10% for gamma-ray energies below 
100 MeV, 5% at 562 MeV, and 20% above 10 GeV. We used a linear interpolation for energies in 
between. 

The foreground templates depend mainly on the CR spectral shapes. We tuned the CR spectra 
of electrons and protons and other nuclei to best describe the gamma-ray sky using the public 
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latitude [deg.] 


(C) 
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Figure 2: (a) The results from a template fit without the 1/E^-^ template to the gamma-ray spectrum in the GP 
for longitudes centered around 35°. The numbers in brackets in the legend indicate the normalization with 
respect to the Dragon prediction, (b) The high energy tail is well described after adding the 1/E^-^ template 
to the fit. The strength of the 1/E^-^ template is indicated in brackets in units of the black solid curve in 
Fig. Id. (c) Latitude distributions for the diffuse component (outer histogram) and the 1/E^-^ emission 
(inner green histogram, normalized to the outer histogram at the center) for —30° < / < 40°. The dashed 
line corresponds to the gas distribution as implemented in the Dragon propagation code. The solid line was 
taken from the molecular column density in Ref. [46]. (d) Latitude distribution of the 1/E^-^ template fiux at 
an energy around 2.4 GeV on a logarithmic scale for —30° < I < 40°. 


propagation code Dragon [12]. The resulting foreground templates are shown in Fig. la-lc. Here 
we superimposed the templates in different directions in a cone of 60° around the GC. All spectra 
were normalized at 3.4 GeV, so only differences in shape, not in flux, are shown. One observes 
that the template does not change with direction, as expected, since the energy losses are small 
for nuclei. The Bremsstrahlung differences are caused by the varying magnetic field in the Galaxy. 
For IC the differences originate form the difference in the spectra of the ISRF, which is donfinated 
by the Cosmic Microwave Background (CMB) in the halo, but has contributions from infrared and 
stellar light in the GP. The template differences in different directions were taken into account in 
the template fit. Towards the GC the isotropic contribution is small and since its normalization 
is fixed mainly by regions outside our region-of-interest and the normalization is the same for all 
subcones, its normalization is fixed in our fit. Varying the isotropic component by ±10% from its 
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Figure 3: (a) Flux from the 1/E^-^ template in the plane. The contour corresponds to the morphology 
found by the authors of the discovery paper [47]. The units correspond to the bubble template unit in Fig. 
Id. (b) Longitude distribution for the flux from the 1 /E^'^ template in the GP(|Z7| < 0.5°). 


nominal FST value for the diffuse class of the P7REP_SOURCE_V15 data worsened the fit. 

Leaving n 3 and n 2 to be both free parameters in the fit leads to unstable fits because of the 
strong correlation of the BR and IC contributions at low energies. Therefore we take the ratio 
n 3 /n 2 of the two templates from the Dragon program for each sky direction. Varying this ratio by 
up to 30% in each direction did not significantly change the results, nor the lower limit of 10% for 
IC or BR of the Dragon flux. A lower limit was imposed to stay away from unphysical fits without 
IC or BR contributions. Einally, there are only 3 free parameters for each subcone, namely ni, n 2 
and 

The energy spectrum in the Eermi bubbles is significant harder [18] and is well described 
by a \/E^-^ proton injection spectrum without any break, as shown in Eig. Id. The depletion 
below 1 GeV originates from the kinematics of production. An example of a template fit with 
foregrounds only is shown in Eig. 2a for a single subcone. One observes that the data below 10 
GeV are well described by the sum of the foregrounds. However, above 10 GeV the data show an 
excess. This is the observation of the l/E^-^ template in the GP. Including this template leads to 
the fit shown in Eig. 2b, which describes well the high energy excess. 

The template fit discussed above was repeated for all subcones. Surprisingly, with only 3 
free parameters for each direction (^ 21 ,^ 22 ,^ 24 ) we can describe the whole gamma-ray sky with an 
excellent except for some parts of the GP. Here we needed to adjust the template towards 
regions with a high column density of MCs, as will be discussed in the last section. The correlation 
with the MCs is demonstrated in Eig. 2c: the distribution in latitude of the l/E^-^ template is clearly 
narrower than the one of the template and corresponds to the latitude distribution from the 
column density of MCs from Ref. [46], shown by the inner solid line. The diffuse component has a 
wider distribution, as shown by the outer dotted line, which was taken from the gas distribution, as 
implemented in the Dragon code. To get agreement between the width of the latitude distribution 
and the gas distribution it was important to take the \/E^-^ template into account and to include 
the breaks in the injection spectra (see Sect. 6 ), since otherwise no good fit is obtained and the 
latitude distribution is far too wide. 
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Figure 4: Sketch of the Milky Way with the central bar and spiral arms. The picture was adapted from [48]). 
The "slim (yellow) ellipse" in the center indicates the angle and length of the region of the 1/E^-^ template 
obtained from the endpoints of the bubble-like emission in Fig. 3b. 


The flux from the template is a factor 40 stronger in the GP than in the halo, as 

demonstrated by the narrow peak in Fig. 2d. In the GP the intensity ratio R of the SRCs and 
production from diffuse Galactic CRs was predicted [16]: 


R(E) = 0.07 


^nSCr^ 


10 ^ yrj \l GeV 


0.6 


(2.3) 


where the first bracket takes into account the difference in gas densities in SNRs and the GP, 
the second bracket the limited confinement time Tp of SCRs in SNRs and the third bracket the 
difference in energy dependence between the bubble template and the template. From Fig. 
2b one observes that this ratio becomes one for energies around 20 GeV, which requires that the 
product of the first two brackets is of the order of a few, which is the right order of magnitude, 
given that the gas density ratio itself is already of this order of magnitude. 


3. Morphology of the Fermi Bubbles 

The fitted values of ^ 4 , i.e. the flux from the \/E^-^ contribution, in all sky directions are 
shown in Fig. 3a, which obviously show a similar morphology as found by the authors of the 
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discovery paper [3] (indicated by the surrounding line), but we find a much richer structure inside 
because of the better spatial resolution with our method. However, the surprise is a strong bubble¬ 
like emission in the GP, shown by the intense (red) bar in the center of Fig. 3a. Its distribution as 
function of longitude is shown in Fig. 3b, which reveals a strong increase at / = +39° followed 
by a sharp decrease at / = — 31°. These angles are close to the endpoints of the bar, as pictured 
in Fig. 4. Assuming these to be the endpoints and taking a distance between the GC and the Sun 
of 8.3 ±0.4 kpc [49] the morphology of the bar is completely determined, as is apparent from the 
geometry shown in Fig. 4. Using a conservative error for the directions to the endpoints of ±1° we 
find from a fit that the major bar axis makes an angle of 77.7° ± 2.1° with respect to the GC-Sun 
line and has a half length of 5.9±0.1 kpc. This morphology is shown by the slim (yellow) ellipse in 
Fig. 4. It has a different angle than the bar indicated in the background, which has 44° ± 10° [48], 
but this angle was found from a subset of stars in the first quadrant. We cannot claim that the 
complete contribution originates from the yellow bar region, since spiral arms before the 

bar may contribute as well. Evidence for contributions from spiral arms comes from the fact that 
we see the bubble-like emission also at longitudes around -50° as a separate peak in Fig. 3b. This 
is exactly the region of the tangent point of the Scutum-Centaurus spiral arm, as shown in Fig. 4. 

We interpret the combined findings as follows: SNRs accelerate CRs with a spectral index 
of 2.1 (on average) by diffusive shock wave acceleration, which in turn produce gamma-rays with 
such a hard spectrum as long as they are connected to the sources, the SCRs. Towards the GC the 
global thermal and CR pressure can be high enough to blow a small fraction of the gas into the halo, 
either from the GP [50,51] or from the inner GC [52,53]. This hot gas in the halo is inferred from 
the ROSAT X-ray data [54], which was interpreted as advection of gas from the GP [55]. In such 
an advective environment the CRs can be trapped inside the plasma without an energy-dependent 
escape, so the SCRs will still correspond to the l/E^-^ spectrum at high latitudes or the CRs are 
accelerated in the shocked gas from the outflow [5]. The spectral index of the locally observed 
proton spectrum is around 2.8 and the softening from 2.1 to 2.8 is attributed to diffusion, simply 
because high energy protons escape faster from the sources and the Galaxy [2]. 

4. Energetics 

The luminosity of the 1/E^-^ emission in the halo between 1 and 100 GeV for 10° < \b\ < 55° 
and |/| < 30° corresponds to 5.6±0.3{fit)±0.9{sys) • 10^^ erg/s. Here we followed the calculation 
in the discovery paper by [3], who found 4 • 10^^ erg/s for the halo bubbles without giving an error. 
The first error originates from the fit, while the second error originates from the error in the spec¬ 
tral index, which is the dominant systematic uncertainty because the normalization of the bubble 
template is most sensitive to the high energy tail in the data, so the extrapolation to lower energies 
depends on the spectral index. The bubble luminosity increases to 8.7 ±0A{fit) ± 1.5(^y^) • 10^^ 
erg/s, if we decrease the latitude limit to 1.5°, while the 1/E^-^ emission in the GP (\b\ < 1.5°) has a 
luminosity of 6.9 ± 0.6{fit) ±1.1 (sys) • 10^^ erg/s. This can be compared with the hadronic energy 
release from SNRs in the inner Galaxy. We expect for the luminosity in gamma-rays between 1 and 
100 GeV: Ey = Ecr % ^(i-ioo) ^scr/{scr+gcr) erg/s. Here Ecr is the fraction of the SNR me¬ 
chanical power ±0 into hadronic CR energy, Ey is the energy converted into gamma-rays, £(i_ioo) is 
the fraction of gamma-ray energy between 1 and 100 GeV, and ^scr/{scr+gcr) is the fraction of en- 
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Figure 5: Morphology of the 1/E^-^ contribution (histogram) in the GP (lbl< 2° in comparison with the 
morphology of the ^^Al line (bottom graph from the Comptel space mission 119] and the insert from the 
Integral space mission 120]). 


ergy into SCRs. Hillas [17] estimates a SNR rate of 1 SNR/century in the inner region (R <4 kpc) 
corresponding to Eq = 10^^ erg/3A0^ ^ = 3.3 • 10^^ erg/s. Using Ecr £y £i-wo % c /?/(^ c /?+ gc /?)= 0.16 
X 0.1 X 0.3 X 0.3=0.0009 we find Ey = 4.8 • 10^^ erg/s for the gamma-ray energy from SCRs to be 
compared with a summed bubble-like emission in the halo and disc of Ey^^ = 1.6 • 10^^ erg/s. This 
is reasonable agreement given the large uncertainties involved. We obtained ^scr/{scr+gcr) = 0-3 
from the energies in the and ni components in the fit, while the other efficiencies were estimated 
as in Ref. [17]. 


5. Comparison with sky maps 

The radioactive isotope ^^Al is synthesized in stars and expelled into the interstellar media by 
winds or supernovae explosions. It can be observed by its characteristic 1.809 MeV gamma-ray 
line emitted in the decay (lifetime 7.10^ years). Given the small j8 = v/c factor in combination with 
its short lifetime the element is expected to decay close to the source and hence to be a good tracer 
of SNRs. Although other stars do contribute, SNRs are expected to be the dominant source. Since 
the SCRs trace SNRs as well, we should expect a close correlation between the skymaps of the 
1.809 MeV signal of ^^Al and the multi-GeV hard component in the Fermi gamma-ray data. This 
is indeed the case, as can be seen from a comparison of the longitude distributions of the \/E^-^ 
and ^^Al fluxes [19,20] in Fig. 5. Note that the spatial resolution of the Fermi LAT instrument is 
around 0.5°, while for the ^^Al data from Comptel it is about 1.6°. 
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(a) 


(b) 




(c) (d) 




(e) 


(f) 


Figure 6: Left: some examples of the template fit in regions with strong ^^Al production without using a 
break in the injection spectrum. The few GeV excess in the data is clearly visible at longitudes even far away 
from the GC (/ = 80°). A break is needed in all directions with strong ^^Al production, see text. Right: as on 
the left, but the excess disappears for a break in the injection spectra, as indicated in Fig. 7. The higher the 
MCs density, the higher the break. Note that the break in the Pamela data is thought to originate from stellar 
winds and energy losses, so in regions of dense MCs it is not unreasonable to have deficit of low energy 
protons. The numbers in brackets in the legends indicate the normalization of the described contribution 
with respect to the Dragon prediction or the strength of the bubble template in units of the solid curve in Fig. 
Id. 
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(a) 


(b) 


Figure 7: (a): Several proton spectra describing best the gamma-ray data in various regions: black dashed 
line for most of the sky, i.e. outside the regions of MCs. The solid line corresponds to the dashed line after 
applying a solar modulation of 154 MeV, which can be compared with the Pamela data (purple circles). 
The dashed lines with depletion at low rigidities are the spectra with a 7 (blue) and 14 GV (green) break, 
respectively. Breaks between 7 and 14 GV are needed to describe the gamma-ray spectra in MCs (see text), 
(b): the gamma-ray spectrum for various slices of the proton spectrum. Note that the maximum photon 
energy is always below the proton energy. One observes clearly a shift of the maximum flux, if the low 
energy protons are suppressed by a break as indicated in (a). 


6. The GC GeV Excess (GCE) 

Towards the GC the template fit gets considerably worse, since the data show an apparent 
excess at a few GeV, as demonstrated in Fig. 6a. This excess turns out to be present everywhere in 
the GP, where the ^^Al production is strong. A few examples are shown in Figs. 6c and 6e. This 
excess can be removed by a break in the injection spectra, leading to a shape of the proton spectra 
in Fig. 7a. The curve at the left side shows the interstellar spectrum best describing the gamma-ray 
sky outside the GC and bubble regions, i.e. for |/| > 40°. This tuning is called "gamma-best". After 
solar winds are taken into account the "gamma-best" proton spectrum resembles the observed local 
proton spectra, as shown by the Pamala data points. Energy losses and stellar winds deplete the 
observed spectra at low energies. It is not unreasonable to expect that these effects deplete the 
lower energy part of the proton spectra even more in the environment of dense MCs. Such a deficit 
can be parametrized by a break in the injection spectrum, as shown by the dashed lines in Fig. 7a 
for a break at 7 and 14 GV, respectively. Such breaks immediately remove the apparent gamma-ray 
excess, as shown on the right hand side of Fig. 6. This can be understood from Fig. 7b, which 
shows how a deficit of low energy protons in the parent proton population shifts the maximum 
in the gamma-ray spectrum from below 1 GeV to a few GeV. A break at 7 GV is needed for a 
good fit towards the tangent point of the spiral arm at a longitude of -50°, but in directions with 
a higher column density of MCs (defined by a higher flux of ^^Al production) the needed break 
shifts to slightly higher rigidities: 8 GV for the "Cygnus" region at 77.5° < / < 82.5°, 9 GV for 
the direction of the bar (—31° < / < 39°) and 14 GV towards the Central Molecular Zone (CMZ) 
at — 1° < / < 1.5° [56,57]. The density in the CMZ is several orders of magnitude higher than the 
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longitude [deg.] 



longitude [deg.] 


(C) 


(d) 


Figure 8: (a) Differences between the fit and the data (for energies around 2.4), if the templates are based 
on the injection spectra from the locally observed CRs. The colour coding indicates the relative difference 
between data and the best fit value for values between 0 and 50%. One observes a clear excess in the bubble 
region of about 30% and the excess near the GP varies between 60 and 120%. (b) As in (a), but instead of 
the relative difference the absolute difference is plotted in units of 10“^ GeV cm“^ s“^sr“^. Because of the 
high fiux in the GP the large relative error causes a large absolute excess (varying from 20 at 111=60° to 140 
at the GC in units of 10“^ GeV cm“^ s“^sr“^) and the excess at higher latitudes is just an artefact of the 
1.0.s crossing the GP with this large excess combined with a higher absolute fiux towards the GC. (c) As in 
(a) but after optimizing the injection spectra of protons and including the template, which removes 

the excess in all regions, especially the strong excess in the GP, which is largely responsible for the excess 
seen in (b). (d) The absolute excess disappears after the excess in the GP disappears with the tuning in (c). 


averaged gas density. 

How does our analysis compare with other analyses discussing the GCE? All previous analyses 
mask the region in the GP, since the diffuse gamma-ray emissivity is here poorly known. However, 
it is a misconception, if one believes not to be sensitive to the GP in this case, since the integral 
along the l.o.s. of the gamma-rays always crosses the GP and any excess in the GP will show up 
as an excess in whatever direction. Since the column density in the GP increases with decreasing 
latitude, the effect becomes strongest at low latitudes. This is demonstrated in Fig. 8 a, which 
shows the relative residuals, defined as (data-model)/data. Here the model is the predicted flux 
for templates from nuclei injection indices without a break. One observes about 30% deviations 
in the regions of the Fermi Bubbles (even at the low energy of 2.4 GeV, for which the residuals 
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are shown), but the deviations in the GP are much larger and vary between 60 and 120%. If 
we plot these deviations as an absolute excess, we observe the famous excess in the l,b-plane 
towards the GC (see Fig. 8b). It is largest in the GP because of the large flux there combined 
with the large relative deviation (60-120%), but it is visible up to latitudes of 20°. The latitude 
dependence is easily explained by the increasing column density of the GP to the integral along 
the I . O . S ., as discussed above. But why the excess is strongest at small longitudes? The answer is 
simple: the flux towards the GC is increasing (because of the higher gas density) and the excess is 
a certain percentage, so the absolute value of the excess increases with the absolute flux, i.e. with 
decreasing longitude. If we remove the excesses in the GP by retuning the injection spectrum (via 
the breaks in the regions of strong , as discussed above) and add the I/E^-^ template, the observed 
and calculated fluxes are in excellent agreement in the whole sky (see Fig. 8c) and the "pseudo-GC” 
excess disappears, as shown in Fig. 8d. Note that a good fit requires to take the l/E'^-^ template 
into account at all latitudes, even at zero degree, as is apparent from Fig. 3a. We only optimized 
the dominant hadronic emission in the GP, which was sufficient to eliminate the GCE. This clearly 
proves that the apparent GCE is due to the uncertain emissivity in the GP, not in the GC. 

7. Summary 

In summary, with our template fitting method we confirm the bubbles in the halo, but find in 
addition strong bubble-like emission in the GP. Both, in the halo and the GP, the spectral index of 
the hard gamma-ray component corresponds to a proton injection index of 2.1 without a break , as 
expected for "freshly" accelerated CRs inside the sources (SCRs). The narrow latitude distribution 
of the 1/E^-^ contribution shows that it originates mainly from the region of the molecular clouds, as 
expected for SCRs. In addition, the longitude distribution of the 1/E^-^ contribution coincides with 
the emission of ^^Al, which is presumed to be a tracer of SNRs [21]. The sharp edge of the bubble¬ 
like emission in the GP at longitude 1=-31° in Fig. 3b) coincides with the endpoint of the Galactic 
Bar, while the bump at l=-50° coincides with the tangent point of the nearby Scutum-Centaurus 
arm (Fig. 4). These three observations (agreement with predicted SCR spectrum, correlation with 
SNRs via the correlation with the ^^Al line and correlation with a nearby spiral arm are strong 
arguments for associating the 1/E^-^ contribution with SCRs. 

Why the SCRs and the Fermi Bubbles have the same 1/E^-^ spectrum? This is most easily 
explained, if one assumes the bubbles to be outflows of the star-forming region near the GC, where 
the combined thermal- and CR pressure is high enough to overcome the gravitational pull. If the 
protons are trapped inside the plasma, their spectra are not softened by energy-dependent escape. 
This leads to the same hard 1/E^-^ gamma-ray spectra inside the sources and in the bubbles. Al¬ 
ternatively, the protons may be accelerated in the shock wave of the outflow, thus having the same 
spectrum as the SCRs [5], which are accelerated in shock waves as well. 

Finally, we showed that the famous few GeV excess towards the GCE also occurs in the GP 
(GPE) and is strongly correlated with the ^^Al flux, just as the 1/E^-^ contribution. As discussed 
above, the latter contribution was interpreted as the contribution from "freshly" accelerated CRs 
inside sources (SCRs). In contrast, the few GeV excess was interpreted as the contribution from 
"old" CRs obtaining a deficit at low rigidity in the dense environment of MCs, characterized by 
higher energy losses and stronger stellar winds. Since the SCRs and dense MCs are both traced by 
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production, it comes as no surprise that the few GeV excess (from "old" CRs) and the 1/E^-^ 
contribution (from "fresh" SCRs) have an identical morphology, namely the one from the A1 flux. 
The reason why the GPE shows up as a GCE is simply that the l.o.s towards the GC cross the GR 
If we retune the injection spectra in the GP to correctly describe the gamma-rays for regions with 
a high column density of MCs, the GCE indeed disappears, as expected if the GCE is mainly an 
artefact of the GPE. 
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